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ABSTRACT 

We find the nine bulk-flow and shear moments from the SFI++ sur- 
vey, as well as for subsamples of group and field galaxies. We constrain the 
velocity power spectrum shape parameter F in linear theory using these mo- 
ments. A likelihood function for T was found after marginalizing over the 
power spectrum amplitude crgrij^^ using constraints obtained from compar- 
isons between redshift surveys and peculiar velocity data. We have estimated 
the velocity noise cr* from the data since without it our results may be bi- 
ased. We also performed a statistical analysis of the difference between the 
field and group catalogues and found that the results from each reflect the 
same underlying large scale flows. We found that we can constrain the power 
spectrum shape parameter to be F = O.lSlg Jg for the groups catalogue and 
F = 0.09^0 04 for the field galaxy catalogue in fair agreement with the value 
from WMAP. 

Subject headings: cosmology: distance scales - cosmology: large scale structure of the universe - cosmology: obser- 
vation - cosmology: theory - galaxies: kinematics and dynamics - galaxies: statistics 



1 INTRODUCTION 

Structure formation is assumed to arise from small Gaussian initial fluctuations (Bardeen et al. 1986; Eisenstein 
& Hu 1998) amplified by gravitational instability into the large scale structure we observe in current surveys. This 
scenario has been supported by observations of the bispectrum of galaxy redshift surveys (Scoccimarro et.al. 2001; 
Feldman et.al. 2001; Verde et. al. 2002) as well as baryonic acoustic oscillations (Cole et al. 2005; Eisenstein et.al. 
2005). Since Gaussian random fields are characterized entirely by their power spectrum, on scales where strongly 
non-linear effects are negligible, the power spectrum provides a glimpse of the primordial fluctuations. On scales 
<^ 100/i~^Mpc (where h is the Hubble constant in units of 100 km s~^ Mpc~^) the peculiar velocity field of galaxies 
and clusters is the one of the most promising and reliable probes of the power spectrum (Strauss & Willick 1995). 
Studies of bulk flows of various peculiar velocity surveys (Pike & Hudson 2005; Park & Park 2006; Sarkar, Feldman 
& Watkins 2007) showed the consistency of bulk flow vectors across independent surveys. In a recent paper (Watkins 
& Feldman 2007) we have shown that there is also a consistency between independent proper distance surveys using 
both bulk flow and shear moments. 

Recently there are more and more observations of the peculiar velocity field. The surveys are deeper, denser, 
and more reliable. As surveys have gotten larger, our understanding of the distance indicators needed to extract 
the peculiar velocities has also improved. We are getting better at correcting for the various Malmquist biases and 
other systematic errors, and are able to extract more and better information from surveys (Feldman et al, 2003b; 
Radburn-Smith et.al 2004; Pike & Hudson 2005; Sarkar, Feldman & Watkins 2007; Watkins & Feldman 2007). In 
this paper we focus on the SFI-f-f surveys (Masters et.al. 2006; Springob et.al. 2007), the largest reliable survey of 
peculiar velocities of the local Universe assembled to date. 
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Peculiar velocities of individual galaxies or groups are dominated by errors that are proportional to the distances 
to those objects. In addition, individual galaxy or group velocities contain contributions from small scales and axe 
not described well by linear theory. Previous work has shown that the analysis of individual galaxy velocities using 
linear theory can lead to biased results (Watkins et al. 2002; Feldman et.al. 2003a). However, large-scale moments 
formed from weighted linear combinations of the individual velocities are expected to be better described by linear 
theory and less susceptible to bias. Moreover, since large scale moments have much smaller uncertainties than 
individual velocities, these moments concentrate the information contained in the survey. In this paper we will work 
with the nine lowest-order moments corresponding to the three bulk flow and six shear moments, thus keeping much 
of the information contained in the survey while avoiding possible biases due to small scale flows. 



2 ANALYSIS 

The velocity field can be expressed in a Taylor series 

Vi(r) = Ui -l-pi^r^- -I- . . . (1) 

where u is the bulk flow vector and pij is the shear tensor. The surveys we consider are a set of A'^ galaxies (or 
groups) each with a position vector r„ and a line-of-sight velocity Sn = v-f„ with measurement error it„. Following 
Kaiser (1991), we create a 9-component vector gi{r) — {fx^fy^f^^rfxfx, rfxry,rfxfz,rfyfy,rfyfz,rfzfz) to model 
the bulk-flow and shear compontents of the velocity field. What follows is similar to the development of Watkins & 
Feldman (2007); see that paper for details. 

Individual survey velocities can be modeled as 

Sn = Opgp(r„) + 5n, (2) 

where Op are the moments and 5n are Gaussian distributed random variable with zero mean that accounts for 
deviations from our model. We can write the variance of 6n as cr^ -|- cr^ , where an is the measurement error and ct* 
represents the small-scale higher order linear and non-linear contributions that were not taken into account by our 
expansion. 

With these assumptions, the maximum likelihood estimate (MLE) for the moments is given by 

JV 

ap = ^ Wp (3) 

m — 1 

with weights given by (Kaiser 1991; Watkins & Feldman 2007) 

Afmom , , 

».,n= (4) 

(J 1-1 ~1~ (J J, 



an + at 

9=1 



where 



A _ gp(rm)g9(rm) . s 

Here iVmom = 3 for the bulk flow only and iVmom = 9 for the bulk flow and and shear moments. 
The covariance matrix of the moments is given by 

(Opftc,) = {C^Wp,nSn)C^'Wq,mSm)) = ^ Wp,„Wq,m (^n'S'm) ■ (6) 

n m n.m 

where Eq. (2) can be used to write the covariance matrix for the individual mcEisured velocities (SnSm) in terms of 
the velocity field v(r) as 

(S'nSm) = (f„ • V(r„) Tm-^{TCm)) +5nni{o-l + an)- (7) 

In linear theory the first term can be expressed as an integral over the density power spectrum i-'(fc), 
(f„ • v(r„) im ■ v(r„)) = t^htl f dk Wmn{k) P{k) (8) 



27r2 

where f{0,m) ~ f^m^ for the scales we are probing and Wmn{k) is the angle averaged window function 
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= / dk Wl,{k)P{k), (11) 



Wmn{k) = j (fn-k • k) exp (ik • (r„ - r^)) . (9) 

Plugging Eq. (7) into Eq. (6) and using equation (8), the covariance matrix of the moments reduces to two 

terms, 

= 49'+-Rp.'- (10) 

The first term is given as an integral over the matter fluctuation power spectrum, P{k), 
Jo 

where the angle-averaged tensor window function is 

y^pqik) = ^Wp,nW^,mWmnik) (12) 

The second term, called the "noise" term, is given by 

Rpq = ^ Wp,nWq,n {<^n + <^l) (13) 
n 

For the case a = b, Eq. (12) gives the angle-averaged window function for the moment Op. This window function 

tells us which scales contribute to the value of the moment. For the MLE weights, the bulk flow and shear window 
functions are determined by the geometry of the survey and the velocity measurement errors. 

Given a peculiar velocity survey and the values of its N^nom moments, we can write the likelihood of a theoretical 
model used to calculate the covariance matrix as 

^ = J^r72 (- ■ (14) 

We use this equation in order to place constraints on the parameters of cosmological models, in particular F, the 
parameter that determines the shape of the power spectrum as will be discussed below. 

In general, the bulk flow probes scales larger than the diameter of the survey while the shear component probes 
scales similar to the survey's radius. In this study, we will use a prior constraint on ugO^® fixes the power spectrum 
amplitude on scales quite a bit smaller than that of the survey. Given that the shape parameter F controls the 
relative distribution of power between large and small scales (smaller F corresponds to relatively more power on 
large scales), increasing F generally results in smaller values for the elements of the covariance matrix. Surveys with 
large values of bulk flow and shear moments generally favor smaller values of F (see Eq. 14). 

How a survey samples the power spectrum is a function of the distribution of the survey objects as well as 
their measurement errors; both of these properties are reflected in the survey's window function. The values of the 
moments are not strictly comparable between surveys since each survey has a unique window function (Watkins 
& Feldman 1995; Sarkar, Feldman & Watkins 2007; Watkins & Feldman 2007). However, given the surveys we are 
considering here, the values of the moments should be highly correlated. To quantify the agreement, we use the 
covariance matrix for the difference between the moments for the two surveys, 

iXpq \V^P }y^Q )/ — ^pq ~r J^pq ^pq ^qp : \ / 

where the cross-terms are given by 

Rpf = dk (W^'')lq{k)P{kl (16) 



r 

i 

(^^)-i(^^)-/ ((ai?)2 + (a?)2j((lS)^ + {am Jl^(^^-^^^-^) '^^P [* • - '■^)] 



and the squared tensor window function is 

(W'^'fpqik) ~ 



(17) 



Here we assume that the nonlinear contributions to the velocities of the galaxies represented by cr, in the two 
surveys are uncorrelated. This is not likely to be true in reality, since galaxies in the same local neighborhood are 
affected by the same small-scale flows. However, this will always cause us to underestimate the expected amount of 
correlation between surveys. Thus our results on how well two surveys agree should be considered as upper bounds. 
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Given that the covariance matrix is a convolution of the power spectrum and the window function, (Eq. 11) times 
f^i^m), we follow Watkins & Feldman (2007) and adopt the Eisenstein & Hu (1998) parametrizing the transfer 
function with the parameter F. 



3 SURVEYS 

The surveys wo apply our formalisms to are the SFI++ (Masters et.al. 2006) catalogue, as well as subsets of groups 
and field gala:x;ies. The entire catalogue consists of ~ 5000 spiral galaxies that have I-band Tully-Fisher distances 
and velocities (Tully & Fisher 1977). The catalogues represent a homogeneously derived set of positions and peculiar 
velocities that arc corrected for both the homogeneous and inhomogenoous Malmquist biases (Springob et.al. 2007) 
using the 2MASS redshift survey (Huchra et.al. 2006; Masters et.al. 2007) for the large-scale structure estimate. 
The final catalogues we use have 2713 field galaxies with an effective depth of ~75ft~^Mpc and 736 groups with 
an effective depth of '^45/i^^Mpc. These samples are based primarily on the all-sky Spiral Field I band (SFI) and 
the Spiral Cluster I band (SCI) samples compiled in the 1990's, (see (Haynes et al. 1999a; Haynes et al. 1999b) 
and references therein) as well as some data from other samples. The sample covers most of the sky with galactic 
latitude b > |10°| with some deficiency of galaxies iu the declination range of [—17.5°, —2.5°]. The rotation widths 
for the SFI-I— I- sample come from both the 21 cm line global profile widths (~60%) and optical rotation curves 
(~40%). 



4 RESULTS 

In Table 1-a we show the MLE for Op for the field galaxies (SFI-I— |-f), groups (SFI-I— |-g), and complete (SFI-I— |-) 
catalogues for the nine bulk flow and shear moments. The uncertainties given are the diagonal terms in the "noise" 
term of the covariance matrix (see Eq. (13)). The covariance between the different term, i.e. the off-diagonal terms 
in the covariance matrix, are at most 10% of the diagonal terms. The moments and the value for a* given for 
each survey are found using an iterative procedure which maximizes the likelihood; see Watkins & Feldman (2007) 
for details. The likelihood distribution for ct, is quite broad, so that its uncertainty is on order 50% of its value. 
However, we have found that our results are fairly insensitive to the precise value used in the analysis. 



Table 1-a Bulk flow and shear magnitude and uncertainties 

Galactic coordinates in the CMB rest frame 



name 

Components 


X 


Bulk flow 
km 

y 


z 


XX 


xy 


km i 

xz 


Shear 
3-iMpc-i 

yy 


yz 


zz 


km s~^ 


SFI-H-Hf 


49± 54 


-260±53 


100±45 


3.2±1.7 


-3.4±2.4 


2.2±1.8 


3.7±1.6 


-2.2±1.9 


2.7±1.4 


720 


SFI-H-Hg 


145±66 


-233±69 


71±47 


4.3±2.6 


-3.1±4.2 


7.1±2.9 


-0.43±2.8 


-5.7±3.2 


-0.73±1.9 


545 


SFI-h-h 


83±42 


-266±42 


73±33 


3.3±1.6 


-3.1±2.1 


3.5±1.6 


2.8±1.5 


-3.0±1.7 


1.8±1.4 


652 



Table 1-b is the same as Table 1-a with except that the moment values are calculated assuming a model that 
includes only the three bulk flow moments. We also show the value for a* calculated for three degrees of freedom 
for each survey. As expected, the a* values found using the bulk flow model are larger than those found using the 
bulk flow plus shear model, indicating that at least some of the variation about the bulk flow model is accounted 
for by shear. However, overall the values of a* that we have obtained here are somewhat larger than those found 
for other surveys (see Watkins & Feldman 2007), including the original SFI survey (Haynes et al. 1999a; Haynes 
et al. 1999b) whose estimate for a, was only 413 km/s. This is at least partly due to the fact that the surveys 
include many objects with very large measurement errors. Large measurement errors tend to broaden the likelihood 
distribution of cr« and lead to larger estimates. This effect is seen in the fact that SFI-I— ho, which combines objects 
into groups and hence has smaller errors than SFI-I— fp, also has a somewhat smaller estimate for u*. Of perhaps 
more concern is the fact that the SFI-|--|- survey contains a significant number of galaxies with anomalously large 
velocities that cannot be accounted for by their measurement errors. The presence of these galaxies in the survey 
also results in a larger cr*. We will be exploring these issues further in a future paper. 
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Table 1-b Bulk flow magnitude and uncertainties 



name 

Components 


X 


Bulk flow 
km 

y 


z 


(T* 

km 


SFI++P- 


70 ± 53 


-281 ± 52 


62 ± 44 


735 


SFI++G 


91 ± 63 


-209 ± 66 


39 ± 45 


554 


SFI++ 


72 ± 41 


-258 ± 41 


50 ± 32 


662 



Eq. (14) gives the likelihoods for the parameters of the model for each survey. The theoretical covariance matrix 
for the velocity moments is completely specified by two parameters; the shape parameter F, and the amplitude 
(Tsf^m^- The amplitude parameter has been strongly constrained by comparisons of the velocity field obtained from 
peculiar velocity data and density fields obtained from redshift surveys, (for a review of these constraints see 
Pike and Hudson 2005). We consider the constraint crgOjJ;^ = 0.45 ± 0.05, which corresponds to that obtained by 
Zaroubi et al. (2002), as being most representative of constraints obtained using this method, and we will adopt 
this constraint as a prior in our calculation of the likelihood function for F. An estimate of as^^ = 0.52 ± 0.06 
obtained directly from the SFI-I— I- survey alone (Masters et.al. 2006) is consistent with this prior. (An interesting 
discussion of ways to estimate erg from peculiar velocity surveys is given by Abate et. al. 2008, sec also Watkins et 
al. (2002) and Feldman et al. 2003a. ) This prior is slightly higher than the value obtained by WMAP three year 
data (Spergel et al. 2007) but is consistent with it to better than 2— a. By marginalizing over the prior constraint 
on fTsfi^^ we are able to calculate the likelihood function for the single parameter F which determines the shape of 
the power spectrum (for more details see Watkins & Feldman 2007). 

In Fig. 1 we plot the likelihood functions obtained from each of the surveys. These likelihood functions are 
asymmetric and have nongaussian tails. In Table 2 we give the maximum likelihood values of F for each survey 
together with the region around the maximum that contains 68% of the probability under the curve. We also list 
the at the maximum likelihood, for nine and three degrees of freedom, where 

= ^ttpilpq ttq. (18) 
V,<1 

These values show that the peculiar velocity data is consistent with the power spectrum model we are considering. 
As can be seen in Table 2 and in Fig. 1, the maximum likelihood value for F are consistent for the three and 
nine DoF, but given the information contained in the six additional inomciits, the nine DoF likelihood is narrower, 
providing tighter constrains on F. We would like to stress that the likelihood functions calculated here take into 
account the entire covariance, including the off diagonal terms. 



Table 2. The mciximum likelihood value of F 
and the x^ for three and nine degrees of 
freedom for each of the surveys. The last 



entry is the results for the composite survey. 

3 DoF 9 DoF 


Survey 


F 




F 




SFI-h+F 


U.U»_Q Q4 


3.13 


U.U»_Q Q4 


8.63 


SFI-h+G 


17+0-26 
-0.10 


3.06 




9.47 


SFI-+ 




2.97 


0.09+:;:;;^ 


8. 81 
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SFI+ + 



SFI++^ 



SFI++, 




o 
o 



o 

o 



Figure 1. The T likelihood function for the SFI++ catalogs. The top panels are the maximum likelihood for the bulk-flow 

moments only, the bottom ones arc the likelihoods including all nine moments (bulk- flow and shear). The left panels arc for 
the complete catalog, the middle ones are the likelihoods for the group catalogue, the field galaxies likelihoods are on the 
right. The inclusion of the shear moments provides tighter constraints on F. 



Although the group and field galaxy surveys yield consistent results for the value of F it does not necessarily 
mean consistency in the actual values of their moments. In order to check for more detailed compatibility between 
the surveys we consider the question of whether the differences between the values of the individual moments of 
the two surveys, — a.p , (for field and group catalogues, respectively) are consistent with that predicted by the 
theoretical models, «.e. are the measurement errors, the velocity noise, and the differences in how each survey probes 
the power spectrum large enough to explain the differences in the moments. To do that we again use a analysis; 
calculating the covariance matrix Rpq'^ of the difference (Eq. 15) we form 



X 



El F G\f r,F-G\-l I F G\ 



(19) 



The calculated in this way does not depend very strongly on F in the region of interest. For simplicity, then, 
we report x^ values calculated for the single value of F = 0.14. Other values of F give similar results. Table 3 gives 
the results of this analysis, which shows good consistency between the catalogues for the favored range of F values. 
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Thus the velocity moments of the surveys that we consider agree not only in magnitude, but also in value, inasmuch 
as they are expected to given measurement errors, velocity noise, and differences in survey volumes. 



Table 3: for tlic^ difforc^iicos botwcx^ii tlio sur\'CA'S. 







x' 


Surveys 


3 DoF 


9 DoF 


SFI++P. - SFI++G 


0.647 


5.59 



Since the surveys are consistent with each other, it seems reasonable to combine them into a composite survey 
which can then be used to obtain the strongest possible constraint on F. Since the different values of a* for the 
various surveys reflect differences in the populations and distance measures between the surveys, we assign each 
galaxy in the composite survey the value of cr* of its parent survey. In Figure 1, we also show the likelihood function 
for r resulting from the composite survey, with the maximum likelihood value being 0.09 ± 0.04 . In Tables l-a,b 
we give the maximum likelihood values for the bulk flow and shear moments for the composite survey. In Table 2 
we present the maximum likelihood value of V for the composite survey together with its associated x^- 

In this study we have used a two parameter model of the power spectrum that is strictly valid only for the 
zero-baryon case, where theoretically V = flrnh. By interpreting our results for F it is possible to include the effects 
of baryons to a first approximation. Sugiyama (1995) has determined that F scales with baryonic density fij as 

F = Qmh exp [-Qb (l + V2h/0,m)] ■ (20) 

The parameters in this formula are tightly constrained by microwave background studies (Spergel et al. 2007); 

specifically, h = 0.732^° °^2. = 0.241 ± 0.034, and Qth^ = Om23ta'ooml, which can be combined to give fib = 
0.0416±0.0049. Plugging these values into Eq. (20) and propagating uncertainties gives the result F = 0.137±0.025, 
which is consistent with our results. While this model is not as accurate as that of Eisenstein & Hu (1998), the 
latter introduces complications in interpretation due to its F having k dependence. We have found that the use of 
the more complicated model in our analysis did not change our results significantly given the precision that can be 
achieved using the available data. 

The composite SFI++ catalogue gives an estimate of the mean bulk flow to have a magnitude of 288 km s~^ ± 
71 km s~^ toward I = 285° ±14° and b = 11° ±10° where I and b are the galactic longitude and latitude respectively. 
This result agrees well with other estimates from the SFI catalogue, such as the one obtained by (Giovanelli et al. 
1998) to have a magnitude of 200 km s"^ ± 65 km s"^ toward / = 295° ± 20° and 6 = 25° ± 20°. It also agrees 
with the results in (Sarkar, Feldman & Watkins 2007) of 330 km s~^ ± 101 km s"^ toward I = 234° ± 11° and 
6 = 12° ± 9°. 

The SFI++ catalogues discussed in this paper are highly consistent with each other, and provide a strong 
constraint of the power spectrum shape F. Since the groups catalogue is shallower, it appears that the flows 
reflected by the bulk flow and shear moments are smaller than that for the field galaxy catalogue. This results in a 
flatter power spectrum and thus larger MLE for F. This seems to suggest that we have not reached convergence of 
the velocity field with current generation of surveys. However, this study supports and strengthens the notion that 
velocity flelds have come of age and provide an excellent dynamic probe of the large scale structure of the Universe. 
The results from this and similar analyses show remarkable agreement between independent catalogues that survey 
various morphologies, selection criteria, geometries, in galaxy clusters, groups of galaxies and individual galaxy 
catalogues. Both our surveying abilities and our analyses tools have matured to the point that we can study flows 
on scales unreachable just a few years ago and provide uniquely dynamical probe to the study of the large-scale 
structure of the Universe. 

Acknowlegements: HAF has been supported in part by a grant from the Research Corporation and by the 
University of Kansas General Research Fund (KUGRF). 
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